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Abstract 

In this paper we describe a check of the implementation of SANC system generated 
| modules into the framework of the WINHAC Monte Carlo event generator. At this 

stage of work we limit ourselves to inclusion of complete one-loop electroweak cor- 
rections to the charged-current Drell-Yan process. We perform tuned comparisons 
of the results derived with the aid of two codes: 1) the standard SANC integrator 
with YFS-inspired treatment of the ISR QED corrections and 2) the WINHAC gen- 
erator, upgraded with the SANC electroweak modules and downgraded to the O(o) 
QED corrections. The aim of these comparisons is to prove the correctness of im- 
plementation of the SANC electroweak modules into the WINHAC generator. This is 
achieved through the presented tuned comparisons. 
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1 Introduction 



The main aim of this work is to implement in the Monte Carlo (MC) event generator 
WINHAC [H the complete O(a) electroweak (EW) corrections delivered by the SANC system 
in the form of the Standard SANC FORTRAN Modules (SSFM) automatically generated 
by the system and to perform a cross check of this implementation by means of tuned 
comparisons of a few distributions with simple cuts. Here we limit ourselves to the 
charged current Drell-Yan-like single W production and use the setup which is rooted in 
the convention of TeV4LHC WS tuned comparisons working group, see Ref. j^j: 

pp — ► W + + X — > £ + is e + X. (1) 

For the description of WINHAC and SANC we refer the reader to the literature: for WINHAC 
to Q and for SANC to [4] and to 00. 

For the case of the charged current (CC) and neutral current (NC) Drell-Yan (DY) 
processes an extended description of the SANC approach can be found in Refs. 0] and 0], 
corresp ondingly. 

For the final state QED radiative corrections WINHAC has been compared with the 
Monte Carlo generator HORACE, both for the parton-level processes and for proton-proton 
collisions at the LHC. Good agreement of the two programs for several observables has 
been found 0]. The comparisons with generator PHOTOS also show good agreement of the 
two generators for the QED final state radiation (FSR) 0]. 

A similar event generator for the Z boson production, called ZINHAC, is under devel- 
opment now. Krakow group also works on constrained MC algorithms for the QCD ISR 
parton shower that could be applied to Drell-Yan processes, see, e.g. Ref. Q. 



Many results of tuned comparison of SANC with several other programs were presented 



for CC case in Ref. [llj and [2[ and for NC case in [12J, showing very good agreement. 
This ensures us in a high confidence of NLO EW SANC predictions. 

In this paper we limit ourselves to presenting the numerical tests of the implementation 
of SANC EW corrections in generator WINHAC, detailed description of the implementation 
itself will be given elsewhere. The paper is organized as follows. In Section 2 we describe 
the setup of the tuned comparisons between SANC and WINHAC. In Section 3 we present 
the results of these comparisons for the total cross sections and various distributions, first 
at the Born level then for 0(a) EW corrections, and finally for a model of purely weak 
corrections. Finally, Section 4 concludes the paper. 



1 SANC is available from the project homepages at Dubna http:/ /sanc.jinr.ru] and CERN 



http://pcphsanc.cern.ch 
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2 Setup of tuned comparisons of SANC and WINHAC 



We use the input parameter set as in Ref. [2j, see also comments after Eq. (4.4.37): 

G M = 1.16637 x 10~ 5 GeV" 2 , a = 1/137.03599911, a s (M 2 z ) = 0.1176, 

M z = 91.1876 GeV, T z = 2.4924 GeV, 

M w = 80.37399 GeV, T w = 2.0836 GeV, 
M H = 115 GeV, 
m e = 0.51099892 MeV, = 0.105658369 GeV, 

m T = 1.77699 GeV, 

m u = 0.06983 GeV, m c = 1.2 GeV, m t = 174 GeV, 

m d = 0.06984 GeV, m s = 0.15 GeV, m b = 4.6 GeV, 



\v ud \ = 


0.975, 


\v us 


| = 0.222 


\v cd \ = 


0.222, 


\v cs 


= 0.975 


v cb \ = \v ts \ 


= \v ub \ = 


\Vtd\ 


= \v tb \ = 



0. (2) 

However, we present the results both in the a(0) and one-loop parametrization 
schemes. 

To compute the hadronic cross section we also use the MRST2004QED set of parton 



density functions [13j, and take the renormalization scale, /x r , and the QED and QCD 



factorization scales, /xqed and //qcd> to be /i 2 = /!q ED = /Wqcd = ^w- 

We impose only detector acceptance cuts on the leptons transverse momenta and the 
charged lepton pseudorapidity 

p £ T > 20 GeV, \t](£)\ < 2.5, £ = e, fi, (3) 

f T > 20 GeV, (4) 

where f \ is the missing transverse momentum originating from the neutrino. 

To simplify the conditions of this purely technical comparison, we do not impose lepton 
identification requirements, as given in Table 4.4.49 of Ref. 0], so we provide "simplified 
bare" results, i.e. without smearing, recombination and lepton separation cuts. We present 
our results only for three differential distributions and the total cross sections, at LO and 
NLO, and the corresponding relative corrections, 5ew [%] = ^cnlo/^clo — 1, for two 
processes: pp — > W + + X — > £ + u^ + X with I — e, /i at the LHC in two schemes: a(0) 
and Gfj,. Moreover, we present the results for some well-defined model of "purely weak" 
corrections 5 wea k, given in Subsection 13.2.3^ for the same cases as for <5ew- 

In our comparisons we use the following W -boson observables: 

• a\y'- the total inclusive cross section of the IV-boson production. 
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the transverse mass distribution of the lepton lepton-neutrino pair. 



The transverse mass is defined as 

= ^2p e T p^(l - cos<p) , (5) 

where p^> is the transverse momentum of the neutrino, and 4> lv is the angle between 
the charged lepton and the neutrino in the transverse plane. The neutrino transverse 
momentum is identified with the missing transverse momentum, j$ T , in the event. 

: the transverse lepton momentum distribution. 



dp T 
da 

d\m\ 



the lepton pseudorapidity distribution 



% = - In I tan ) , (6) 



2 

where the lepton kinematical variables are defined in the laboratory frame. 



One should emphasize an important difference between the conditions of these com- 
parisons and that of TeV4LHC WS concerning the subtraction of initial quark mass 
singularities. Instead of the commonly adopted MS or DIS subtraction scheme (as, for 
example, in Ref. 0]), we use here an YFS-inspired subtraction method 



^YFS 

da 



iSR {s,m d ,m u ;e) 



d&fsR n (s, m d, rn u ; e) S^(s, m d , m u ; e), 



(7) 



where 



+QI 



with 



7l~* 

~6 

2uj 



Q u in — 

ml 
3 



lne 

+ 1 



7l~ 

~3 



(9) 



being the dimensionless soft -hard photon separator (cu is the photon energy). The Q u , Q d 
are the electric charges of the up-type and down-type quarks in the units of the positron 
charge and m u , m d are their masses, while s is the centre-of-mass energy squared of the 
incoming quarks. 

Simultaneously, we subtract in a gauge-invariant way the contribution of the ISR hard 
photons, derived using the W propagator splitting technique 15j. In this way the initial 
quark mass dependence drops out from the one-loop level observables. 
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In order to define our "weak" corrections, we will need the YFS corrections for the 
'initial-final" interference 



-<2 

7T 



S S 

Q d In — - Q d In hi 

—t —u 



In 



s-hP w f + MlT\ 



(10) 



L w L w 



+Qd 
Qu 

and for the "final state radiation" 
5jH(s,mi;e) = 



l ln i_( ln J. + 1 , +Ll2 M. ( 



1, s / s 

-In In hi 

2 —u \ —u 



Li 2 1 + - 

u 



7T 

If- 2 



Ot \ ( S \ , 3 , S 7T^ 

[hx—y-2)hxe+-rhi—y- — 



where s, t, u are the standard Mandelstam variables for the parton-level process and mi 
is the charged lepton mass. 



W + -> ^i*v : cc-scheme 




W* -> n*v : G -scheme 
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Figure 1: The Born distributions of Mp from SANC (red diamonds) and WINHAC (solid 
lines) in two schemes and their relative deviations 5 - 
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3 Numerical results 



In this section we present the numerical results of the tuned comparisons between SANC 
and WINHAC, first the the Born level (LO) and then including the 0(a) EW corrections 
(NLO). At the end of this section we compare also the so-called "purely weak" corrections 
which are the difference between the EW corrections and the "QED" corrections defined 
by the terms given in Eqs. (TTOMTTT) plus the corresponding hard-photon contributions. 



3.1 Comparisons at tree level, LO 

We begin with the comparisons at the Born level. In Figs. [TH3] the distributions are 
shown for all three observables under consideration only for fi + final state but in the 
both schemes: a(0) and G M . The lower parts of the figures shows the relative deviation 
A = (W — S) /W between the two calculations (W for WINHAC, S for SANC). 




As seen, the relative deviations lie within the 1 per-mill band, wherever the cross 
section is not very smalH. We do not show the comparisons for electron channel, since at 
tree level the muon mass effects are negligible, and the plots look identical. 

2 On the SANC side we have both a VEGAS [II] based integrator and a FOAM [jj} based event 
generator. In this comparison the integrator has been used. 
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3.2 Comparison at one-loop level, NLO inclusive cross sections 

Turning to the NLO results, we show, first of all, in Table[T]the comparisons of the inclusive 
cross sections (in pb) within the acceptance cuts and the relative radiative correction factor 
(in %), as seen by two calculations (second and third rows). In the first row we show 
SANC results in the conditions of TeV4LHC WS. The numbers agree with those published 
in [2| within statistical errors. 

The Born cross sections from SANC and WINHAC agree well within statistical errors (< 
10~ 4 ). The EW NLO cross sections agree not worse than within a half a per mill or agree 
even within statistical errors in both schemes, both for the electron and muon channels, 
better for the muon channel where we observe the agreement within the statistical errors. 



3.2.1 NLO distributions: electron channel 

We begin the comparisons of the distributions for the electron channel in two schemes for 
our three W observables (M^,p^ and |^|, Figs.HHHl correspondingly) with the "simplified 
bare" cuts. The two upper figures show the quantity 5ew m %> while the two lower figures 
show absolute deviations A = W — S between the two calculations. 

As seen, the 0(a) EW correction <5 EW is quite large (mainly due to the FSR QED 
contribution), it varies by 18% depending on the scheme. It is shifted to the larger 
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LHC, pp 


W + + X 


-> e + u e + X 










a-scheme 


G M -scheme 




LO [pb] 


NLO [pb] 


Sew [%} 


LO [pb] 


NLO [pb] 


Sew 


'%} 


SANC-MS 


5039.19(2) 


5139.33(5) 


1.987(1) 










SANC-YFS 


5039.19(2) 


5137.53(3) 


1.952(1) 


5419.18(2) 


5208.48(3) 


-3.8 


38(1) 


WINHAC 


5039.06(11) 


5138.04(16) 


1.966(3) 


5419.04(12) 


5209.04(12) 


-3.8 


74(3) 






LHC, pp -> 


W + + X 


-> it+v^ + X 










a-scheme 


G M -scheme 




LO [pb] 


NLO [pb] 


Sew [%] 


LO [pb] 


NLO [pb] 


Sew 


'%} 


SANC-MS 


5039.20(2) 


5229.58(6) 


3.778(1) 










SANC-YFS 


5039.20(2) 


5227.73(2) 


3.741(1) 


5419.19(2) 


5305.47(3) 


-2.098(1) 


WINHAC 


5039.03(11) 


5227.87(14) 


3.745(2) 


5419.01(12) 


5305.59(14) 


-2.094(2) 



Table 1: The tuned comparisons of the LO and EW NLO predictions for aw and Sew 
from SANC and WINHAC for the simplified bare cuts. The statistical errors of the Monte 
Carlo integration are given in parentheses. 
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Figure 4: The EW NLO distributions of from SANC (red diamonds) and WINHAC (solid 
lines) for the electron channel in two schemes and their absolute deviations A = W — S. 



negative values in the scheme and more moderate in the a(0) scheme, the reason for 
which the latter was preferred by tuned group of TeV4LHC WS. The absolute deviation 
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for both schemes does not exceed 0.1% in the important regions where the cross section 
is large. For the distributions, it varies within 25% but this is an artificial result of 



W + -» e*v e : a-scheme 



W + e*v e : G LI -scheme 
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Figure 5: The EW NLO distributions of p l T from SANC (red diamonds) and WINHAC (solid 
lines) for the electron channel in two schemes and their absolute deviations A = W — S. 



applying "simplified bare" cuts. The r\i distributions are flat and show little biases of the 
order of a quarter of a per mill. However, most likely VEGAS errors are underestimated 
in the SANC results. 



3.2.2 NLO distributions: muon channel 

We continue the comparisons for muon channels in two schemes for the same three W 
observables (Mp and \rjg\) with the "simplified bare" cuts. The results are presented 
in Figs. [7HH1 respectively. Again, the two upper figures show EW NLO correction 5 EW in 
%, and the two lower figures show absolute deviations W— S between the two calculations. 
Here the absolute deviations in statistically saturated regions do not exceed 0.05% and in 
average is of the order of 0.025%. For the muon channel both calculations are statistical 
consistent and no evident biases are observed. 

It is important to emphasis that biases could be present, in principle, due to finite muon 
mass, which treatment in two calculations is not identical: for the muon channel SANC 
uses fully massive formulae for all contributions while WINHAC uses a mixed approach - 
electroweak virtual and soft real-photon corrections are calculated in the massless fermion 
approximation, while massive fermions are kept in hard real-photon radiation. 
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Figure 7: The EW NLO distributions of from SANC (red diamonds) and WINHAC (solid 
lines) for the muon channel in two schemes and their absolute deviations A = W — S. 
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W + -» |i*v„: a-scheme 



W* -> n'v,,: G M -scheme 



"lo 



: 6 EW [%] 


- WINHAC 




* SANC 
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Figure 8: The EW NLO distributions of p e T from SANC (red diamonds) and WINHAC (solid 
lines) for the muon channel in two schemes and their absolute deviations A = W — S. 
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Figure 9: The EW NLO distribution of \r]i\ from SANC (red diamonds) and WINHAC (solid 
lines) for the muon channel in two schemes and their absolute deviations A = W — S. 
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3.2.3 Weak corrections 



Here we discuss the "purely weak" corrections which are defined as 

^weak 

where 



"softvirt 



eYFS 
"softvirt ' 



rYFS 
"softvirt 



rYFS , rYFS 
°ISR + °Int 



eYFS 
°FSR i 



(12) 
(13) 

with three contributions given by Eqs. (l9lfT0HTT|) . The contribution Sf^ vilt includes the 
1-loop EW corrections plus the real soft-photon correction and is provided by the SANC 
modules. This definition is free of any regularization scales. 

From the Table [2] one sees, that for the electron channel the agreement is very good, 
while for the muon channel we observe the systematic differences of about 0.007%. This 
can be attributed to different treatment of the muon mass in the two programs: SANC uses 
the fully massive formulae while WINHAC uses the massless-lepton approximation for these 
corrections. The "weak" corrections in the a-scheme are quite sizable, ~ 6%, because of 
the light-fermion loop contributions, ~ ln(s/m|), to the W self-energy correction. Such 
contributions drop out in the G^-scheme making the "weak" corrections much smaller, 
~ 0.1%. 



#weak [%, 




LHC, pp -> W + + X -> e+z/ e + X 




a-scheme 


Gfj,- scheme 


SANC 
WINHAC 


5.7223(2) 
5.7220(3) 


-0.1175(2) 
-0.1177(0) 


LHC, pp -»• W+ + X -> + X 




a-scheme 


G^-scheme 


SANC 
WINHAC 


5.7286(2) 
5.7220(2) 


-0.1109(2) 
-0.1177(0) 



Table 2: The tuned comparisons of the "purely weak" corrections <5 wea k from SANC and 
WINHAC for the simplified bare cuts. The statistical errors of the Monte Carlo integration 
are given in parentheses. 



In Figs. [TUHT5l we show the distributions of the "weak" corrections and absolute devi- 
ations between the two calculations. The figures show agreement at the level 0.01%. In 
some cases the biases of the same order are seen. Again, this might be a consequence of 
underestimation of errors by VEGAS. In the muon channel, the observed deviations at 
the level of 0.01% can be attributed again to different treatment of the muon mass in the 
two programs. 
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Figure 10: The "weak" correction distributions of M^f from SANC (red diamonds) and 
WINHAC (solid lines) for the electron channel in two schemes and their absolute deviations 
A = W-S. 
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Figure 11: The "weak" correction distributions of from SANC (red diamonds) and 
WINHAC (solid lines) for the electron channel in two schemes and their absolute deviations 
A = W-S. 
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Figure 12: The "weak" correction distributions of \r)^\ from SANC (red diamonds) and 
WINHAC (solid lines) for the electron channel in two schemes and their absolute deviations 
A = W-S. 
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Figure 13: The "weak" correction distributions of M^f from SANC (red diamonds) and 
WINHAC (solid lines) for the muon channel in two schemes and their absolute deviations 
A = W-S. 
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W* -» ji + v : a-scheme 
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Figure 14: The "weak" correction distributions of from SANC (red diamonds) and 
WINHAC (solid lines) for the muon channel in two schemes and their absolute deviations 
A = W - S. 
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Figure 15: The "weak" correction distributions of |?%| from SANC (red diamonds) and 
WINHAC (solid lines) for the muon channel in two schemes and the absolute deviations 
A = W-S. 
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4 Conclusions 



The main priority of the development of SANC as a HEP tool for the LHC is to create the 
SSFM for the EW corrections at one-loop level to be used in existing MC event generators. 

The goals of this work were: (a) to integrate CC DY SSFM into the Monte Carlo 
event generator WINHAC and (b) to check thoroughly the stability of numbers for simple 
distributions by comparisons of the WINHAC generated results with those provided by the 
recently created SANC CC DY integrator. In this paper we have concentrated on presenting 
the numerical tests of the implementation of the above EW corrections in WINHAC, while 
the details on this implementation will be given elsewhere. 

The main and very important conclusion of this paper is that we have reached the 
agreement between the WINHAC MC event generator and the SANC MC integrator for the 
0(a) EW corrections to the charged-current Drell-Yan process at the sub-per-mill level, 
both for the inclusive cross section and for the main distributions. Thus, our above goals 
have been achieved. 

Another important conclusion is that the MC event generator WINHAC can now be used 
for precision simulations of the charged-current Drell-Yan process at the LHC including 
the 0(a) EW corrections. It can also serve as a benchmark for testing other MC programs 
for this process. 

The next step on this road would be a similar implementation of the SANC modules 
in the neutral-current Drell-Yan MC event generator ZINHAC, being under development 
now. 
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